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Interferometry is an experimental paradigm underlying a range of applications from atomic clocks to gravi- 
tational wave detection and cell-membrane dynamics. The fundamental limit on estimating path differences in 
interferometers is set by the resources available - the volume of space-time, energy, number of particles - through 
the Cramer-Rao bound. Different metrological schemes are judged by the scaling of the estimator's uncertainty 
in the resources consumed. Quantum probes are known to provide enhanced precision in single-parameter es- 
timation. However, multiparameter quantum interferometry remains largely unexplored. Here we show that 
estimating one parameter with quantum-limited precision inexorably leads to a reduced precision of the other. 
Unlike single-parameter estimation, it is impossible to design a strategy saturating the quantum Cramer-Rao 
bound for loss and phase estimation in a single setup simultaneously. We design optimal quantum states achiev- 
ing the best possible simultaneous precisions. Our results reveal general features about concurrently estimating 
hamiltonian and dissipative parameters, and can be applied to sophisticated sensing scenarios such as quantum 
imaging. 



A vast majority of metrological problems, particularly 
those associated with hamiltonian dynamics, can be recast 
into one of phase estimation [1, 2]. In the optical domain, 
high accuracy measurements of the path differences, charac- 
terized by a phase difference tj> between the two modes of the 
apparatus, have found use in many fields [3-6], and has con- 
sistently featured as an important tool in physics, from the 
null result in the search for luminiferous aether [7] to detect- 
ing small variations in the refractive index of biological so- 
lutions [8-11]. Most of these experiments estimate a single 
parameter, the phase difference [12-16]. Quantum states of a 
few photons promise enhanced precision in single parameter 
estimation compared to the best classical strategies [17, 18]. 
In this case, the quantum limit for estimation is always at- 
tainable [13]. This enhanced quantum precision is reflected 
in the improved scaling of the variance of the estimator with 
the number of particles, which is the primary resource in most 
quantum metrological schemes. 

Phase shifts are generated by Hamiltonian evolution of the 
particles. The most general evolution however, comprises 
both Hamiltonian and dissipative dynamics. In an interfer- 
ometric setup, this requires the estimation of phase and loss 
parameters that characterise these two elements of the evo- 
lution. The problem of estimating multiple parameters for 
such evolution is still nascent [19-24]. It clearly has important 
technological implications, not only because all practical sys- 
tems exhibit loss, but also because there are many situations 
in which multiple parameters of several kinds are the objec- 
tive of the experimenter. This is the case, for instance, where 
both dispersion and absorption profiles of a sample are sought 
with high accuracy using a single experimental set up. While 
one can estimate the parameters independently, this demands 
the preparation of different optimal probe states and measure- 
ments for each parameter. This is tantamount to arranging a 
different experimental setup for each parameter and requires 
a full reconfiguration, which may be technically demanding. 
Furthermore sequential estimation is not effective for sensing 
time-varying samples. These limitations are overcome by si- 
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FIG. 1. Schematic of an imaging system : A two-mode n-photon 
quantum probe state | "0o} irradiates the sample and provides the ref- 
erence arm of a two-port interferometer. The detection is chosen so 
to estimate simultaneously the phase shift <j> (Hamiltonian dynamics) 
and the transmission r\ (dissipative dynamics) from the state p(</>, r\) 
resulting from the interaction of the probe with the sample. 



multaneous estimation of both parameters. 

In this Letter, we analyse this measurement strategy, and 
derive fundamental limits for the precision with which phase 
and loss in an interferometer can be estimated in a single ap- 
paratus. We show that it is not possible to attain the individual 
quantum limits of the two parameters at the same time. The 
optimal strategy involves a trade-off between the individual 
quantum precision bounds, revealing how quantum mechan- 
ics forbids their simultaneous saturation. The practical task of 
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accomplishing the optimal precision requires the identifica- 
tion of appropriate quantum probes. To that end, we provide 
the design of optimal two-mode quantum states with a fixed 
number of photons that gets closest to the ultimate quantum 
limit of simultaneous phase and loss estimation, by explicit 
numerical optimization for low photon numbers. For samples 
with high losses, it is natural to use states with a high particle 
flux, and in this limit we present an analytic expression for the 
profile of the optimal quantum probes. 

The archetypal schema for quantum sensing is illustrated 
in Fig. (1). An object, characterized by a set of parameters 
9 = {#„} , is placed in one of the arms of a Mach-Zehnder 
interferometer. This extends the simple case in which a single 
phase shift contains the only relevant information, to a more 
realistic case, in which, for example, both the phase shift and 
the loss are important. The initial probe evolves upon propa- 
gation through the system, acquiring a form that depends on 
the parameter set. At the output, measurements of the probe 
state provide a multivariate probability distribution that cap- 
tures changes in the state due to changes in the system param- 
eter. The optimal probe state is one for which the derivative of 
the probability distribution with respect to these parameters is 
largest. 

For a single parameter, it is possible to identify the optimal 
measurement by analysing the quantum Fisher information. 
This is constructed using the symmetric logarithmic deriva- 
tive (SLD), which quantifies the changes of the state with re- 
spect to the system parameters, and also provides a means to 
identify the optimal probe states and measurement strategies 
that saturate the quantum limit to precision, which is given 
by the quantum Cramer-Rao bound. [13]. In the multipa- 
rameter case, the variation with each parameter can be mea- 
sured by SLDs corresponding to each of the different param- 
eters. The information concerning the parameters 6 that can 
be obtained from the set of any and all possible measurements 
is then encapsulated in the quantum Fisher information ma- 
trix [25]. The covariance matrix, the multivariate version of 
the variance of the estimator, is bounded by the inverse of 
the quantum Fisher information matrix through the quantum 
Cramer-Rao bound [25]. 

The quantum Fisher information matrix involving a pair 
of parameters generating Hamiltonian and dissipative dynam- 
ics is necessarily diagonal, because the tangent spaces corre- 
sponding these two types of dynamics are orthogonal. The 
corresponding covariance matrix for such a pair of parame- 
ters is diagonal when the quantum Cramer-Rao bound is satu- 
rated. The matrix inequality that defines the quantum Cramer- 
Rao bound does not, however, take into account the non- 
commutativity of quantum operators, which may prevent it 
from being saturated (See Supplementary Information). A di- 
agonal covariance matrix implies the statistical independence 
of the parameters. This general feature suggests that Hamil- 
tonian and dissipative parameters may be estimated simulta- 
neously at quantum limit, with their respective uncertainties 
given by the inverse of the diagonal elements, as is the case 
for the quantum Fisher information matrix for phase and loss 



estimation. Yet, it is not so, as we show in the next sec- 
tion. This makes multiparameter quantum metrology nontriv- 
ial, and fundamentally different not only from classical mul- 
tiparameter estimation, but also quantum single parameter es- 
timation, as in both these cases the Cramer-Rao bound can 
always be saturated. 

The best probes states are the ones that maximize the quan- 
tum Fisher information. Since particle number being a central 
resource in quantum metrology, we specify the probe states 
as the most general normalised pure states of a fixed particle 
number, which are of the form 
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(1) 



The covariance matrix for the phase shift 4> (ranging between 
and 2ir) and loss ?/ (ranging between for complete ab- 
sorption and 1 for complete transmission of the probe light) 
associated with the Hamiltonian and dissipative components 
of the states evolution is diagonal when the quantum Cramer- 
Rao bound is saturated. In particular, the quantum Fisher in- 
formation matrix is given by 
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where 1^,1^ are the quantum Fisher information for the 
estimation of phase and loss respectively. Defining s = 
J2t=s x tKt r ) — Y^t=o x tt r ') being the moments of the 
coefficients of Xt, and 6* is a binomial factor [26], 




(3) 



It is easy to conclude the form of the optimal states for the 
estimation of loss and phase independently. In the absence of 
losses, r) = 1, b\ = S sfl , Z^ _1 ^ = 4 (H 2 - H?) , which is the 



variance of X). , maximised for a nOOn state, when X' 
n 2 , as is well-known. In the lossy case, the best states for 
estimation of phase are arrived at by maximizing the general 
form of 2^,0 [15, 27, 28]. On the other hand, the best state 
for estimating the loss is the Fock state \n, 0), in which case 
X nv ~ n [29, 30]. For 77=0, 1 the quantum Fisher information 
for loss X, m diverges. This is also to be expected as in these 
cases all the photons are lost or all transmitted, so that the 
variance in the outcomes of measurements of particle number 
in these cases will be zero. 

The necessary condition for saturating the multiparameter 
quantum Cramer-Rao bound, derived using techniques of lo- 
cal asymptotic normality, is given by the expectation value of 
the commutator of the two SLDs [31]. In our case, this re- 
duces to 

w 

This implies that the optimal measurements necessary to at- 
tain the quantum limits for the two parameters do not com- 
mute, and it is impossible to estimate the phase and loss using 



Tr^L^Z^]] = - f- 



(4) 




FIG. 2. Coefficients of the optimal probe states of the form in Eq. (1) 
for n — 6 obtained by numerical optimisation. The state resembles 
the optimal probes for phase estimation in the presence of losses [15], 
and interpolates towards the Fock state |6, 0) , optimal state for esti- 
mating r\. 



quantum probes with a fixed photon number. The only in- 
stance when the quantum Cramer-Rao bound can be saturated 
is when we learn nothing of the phase <fi. One could still es- 
timate the loss parameter 77 in this situation, using as probe 
the Fock state |n, 0). Evidently, this state and the correspond- 
ing SLD has no sensitivity to the phase in the interferometer. 
More generally, no phase information can be extracted using 
a fixed photon number state as in Eq. ( 1 ) and measurements 
derived from the loss SLD [26]. 

The only option now left is to use the SLD for phase esti- 
mation to estimate the loss parameter. This at least guarantees 
the quantum limit for one of the parameters, while any other 
measurement will be suboptimal for both. Having thus cho- 
sen the SLD for phase to estimate the loss, the classical Fisher 
information is [26] 
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(5) 



This result brings to light another interesting trade-off. It 
shows that the precision of loss estimation can only be en- 
hanced at the cost diminished precision in estimating the 
phase, and vice-versa. In spite of this tradeoff, there is an 
asymmetry between the two parameters. Estimating loss at 
the quantum limit requires a complete sacrifice in estimat- 
ing phase, but phase estimation at the quantum limit leaves us 
with some information about loss. Thus, quantum mechanics 
limits the estimation of Hamiltonian and dissipative parame- 
ters at the quantum limit differently, given by Eqns. (3) and 
(5). 

In addition to the conceptual understanding of quantum 
trade-offs and limitations that this analysis provides, the prac- 
tical challenge in any multiparameter estimation scenario such 
as imaging in Fig. ( 1 ) lies in identifying the quantum probes 
that maximize the precisions of both the parameters simul- 
taneously to the best possible extent. For a single parame- 
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FIG. 3. Comparison of the coefficients of the optimal states from the 
analytical and numerical procedures, for n = 40, 77 = 0.02, w = 
0.9. The red solid line corresponds to Eq. (6). 



ter, maximizing the quantum Fisher information is tantamount 
to minimizing the uncertainty. For multiparameter quantum 
metrology, the sum of the variances of all the parameters is 
bounded from below by the trace of the inverse of the quantum 
Fisher information matrix. In our case of a diagonal Fisher in- 
formation matrix, the optimal probes are then given by the 
minimization of the sum of the reciprocals of and I m 
given by Eqns. (3) and (5). For low particle numbers, this 
minimization can be performed numerically, and the result 
for 71 = 6 is presented in Fig. (2). An alternative criterion of 
choosing the optimal probe states is to maximize the Fisher in- 
formation corresponding to the two parameters. More gener- 
ally, their convex combination which assigns different weights 
to the different parameters can be maximized, and we do so in 
a novel context next. 

Going beyond the current state-of-the-art of few-particle 
quantum probes in quantum metrology [11, 18], we look to the 
future of quantum imaging and sensing involving high-loss 
samples and scenarios such as highly scattering media [32] or 
long-distance quantum interferometers [33]. In these cases, 
high flux probe states are necessary and an explicit numerical 
optimization inadequate. An asymptotic analytic form of the 
coefficients in this case is given by [26] 

1 .9/1 k 
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where 7 a = T",, 1 ^ , and Co is the smallest zero of the first 

' n(l— 77) w ' ^ u 

derivative of the Airy function Ai(y). w is the weight assigned 
to the phase parameter, and in this limit, the precisions of esti- 
mating the phase and loss when using the optimal state scales 
as [26] 
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where (3 = (2w/r](l — ij)(2w — 1)) 2//3 . The linear scaling 
of the quantum Fisher information with the number of pho- 
tons in the probe state for the estimation of phase is the best 



4 



possible in the high loss regime [34]. In Fig. (3), we plot the 
asymptotic profile of the probe states for a high flux, high loss 
scenario, and compare it with explicit numerical computation 
of the coefficients. The similarities between the two profiles 
shows the qualitative and quantitative equivalence between es- 
timating different sets of parameters using the measurements 
derived from the phase SLD . 

We have shown how quantum mechanics limits the simulta- 
neous estimation of multiple parameters at the quantum limit. 
In particular, we have illustrated that the estimation of a pair 
of parameters describing the Hamiltonian and the dissipative 
evolution of a probe state in an interferometric sensor can- 
not give both parameters at the simultaneous quantum limit of 
each. This situation is relevant to a broad range of sensors, and 
opens the door to more complex quantum imaging devices, 
in which multiple parameters related to the object configura- 
tion are sought. Although the saturation of the multiparameter 
quantum Cramer-Rao predicts a diagonal covariance matrix, 
suggesting statistically independent parameters, we show that 
they cannot be estimated concurrently at the quantum limit. 
Furthermore, the best possible strategy imposes a trade-off 
between the attainable precisions of the two parameters. This 
trade-off forces an optimization strategy to determine the form 
of best quantum probes for multiparameter quantum metrol- 
ogy at the ultimate attainable limit. 
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SUPPLEMENTARY INFORMATION 

This supplement discusses in greater detail the technical as- 
pects of the derivation of the results in the main text. 
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Quantum Cramer-Rao bound 

The limit to the ultimate attainable precision of a set of pa- 
rameters 9, their covariance matrix Gov (9), is given by the 
quantum version of the information inequality [25] 



Cov(0) > (X 9 ) 



(8) 



where Xg is the quantum Fisher information matrix associated 
to the quantum state after evolution. This is a matrix inequal- 
ity, and provides as a special case a lower bound on the achiev- 
able error for estimating any single parameter, which given by 
the Cramer-Rao inequality (A6„ v ) 2 > X~^ in terms of the 
quantum Fisher information X vv [13]. For a single parameter 
9 U , the quantum Cramer-Rao bound is always attainable, with 
the optimal measurement being given by the eigenvector of 
the symmetric logarithmic derivative (SLD) defined as 



L v p{9) + p(9)L v = 2 



gg(g) 

d6 v ■ 



(9) 



which also defines the quantum Fisher information matrix as 
T-iiv = 3?Tr [p(9)L l j r L I/ ] . The precision in the estimate of 
the parameters is given by the saturation of the inequality in 
^(AeU) 2 = Tr[Cov(0)] > Tr^Ze)- 1 ] . The classical 
multiparameter information inequality can always be achieved 
asymptotically using unbiased estimators, but in general, the 
quantum matrix equality in Eq. (8) cannot be achieved. This 
is due to the fact that the different SLDs L v , are not gener- 
ally compatible, i.e. they do not satisfy the condition for weak 
commutativity [31] 

Tr\p(e)[L v ,L ll ]]=0. (10) 



Quantum Fisher information matrix 

The orthogonal subspaces which the block diagonal subma- 
trices of 



1=0 



occupy remain orthogonal under infinitesimal translations. 
Therefore, the SLD is similarly block diagonal with (see Sup- 
plementary Information) 



N 



£, = 04 with 



(11) 



1=0 



U v =\du logp^X^I + 2\d v tl>i)(tl>i\ + 2\^)(dM 



Given this block diagonal structure, the elements X^ u of the 
quantum Fisher information matrix are given by 



The first term J M „ [p] = Ya=o Pi d» ^gpid^ log pi is the clas- 
sical Fisher information matrix of the probability distribution 
p. The second term, in some sense the 'purely quantum' con- 
tribution is the weighed quantum Fisher information of the 
constituent states X^ipi^ipi |] = ^Ym=oPi{^ p ^,i)^ where 
we use the shorthand 



Ptiu,k = (d^ k \U k \d v tp k ), 



(13) 



and the operator life = I—\ipk) (ipk | projects onto the subspace 
orthogonal to | ifi^) ■ F° r loss and phase estimation using states 
of the form Eq. (1), the quantum Fisher information matrix 



Ap\ = 



-■r}4> 



x, 



(14) 



is such that X^ = X^ = 0. The diagonals given by Eq. (3). 
Finally, the saturation of the multiparameter quantum Cramer- 
Rao bound, governed by Eq. (10), reduces to 



Tr [p[L n ,L <l> ]}=iJ2Pi{^(d^l\d v ^i}), (15) 
;=o 

which in our case is, in general, nonzero as shown in Eq. (4). 
When estimating both the parameters using as measurement 
the eigenvectors of the SLD for phase La, the precision for 
phase estimation scales at the quantum limit given by 1^. 
That of estimating the loss parameter is governed by the clas- 
sical Fisher information [26] I m = Y^a=oP 1 (®vPl) • 



Fixed-photon-number states 

Consider pure quantum probe states with fixed resources, 
which in our case is the photon number, to estimate simultane- 
ously the loss and phase in the interferometer in Fig. (1). We 
model the loss as a beam splitter of transmissivity 77, resulting 
in a loss mode \1)l- The loss transforms the basis vectors as 



|fc, n — k) <— > y~] 



1=0 



b?\k 



I, n — k) ® where the co- 

1 oic given uy m — y ^ ^ T] k '(1-77)'. Including 
the phase accumulation, the final state is given by 



efficients bf are given by bf = 



n k 

E^ fc0 E 

k=0 1=0 



>bf\k- 



l,n-k)® 



(16) 



with ^fxk = ah- As the information of the photons lost 
cannot be recovered, tracing over the mode L leads to p = 
T1iI=qPi\'4 , i){'4 , 1\i where the normalised states are given 
by 



1 



71—1 
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k=l 



OLk' 



/bfe ik<t, \k-l,n-k). 



(17) 



Owing to the orthogonality of the states (ipi\ipi') — Su>, the 
density matrix can be represented as a direct sum as [27] 



vw = v [p\ + E^v o^x^o ■ (12) 



p 



= (&pi\m^i 



(18) 



1=0 
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Furthermore if the state vectors remain orthogonal un- 
der infinitesimal translations 9 K n- B K + d6 K , that is (i ^ 
j)(d K ipi\ifjj) = 0, and the derivative d K pe and L K decompose 
into diagonal blocks d K pig and L i K supported on the same 
orthogonal subspaces as the pi.g. This decomposition entails 
that eigenvalues of L K are those of the blocks L l K . Each block 
has two non-zero eigenvalues 



A,- K . — 



1 



- i 9 K \ogp l ± J (d K logp,) 2 + 16P i}KK 



with corresponding eigenvectors 

X i , K ^/l\Z\<ip i ) + 2P i \d K ip i 



Pi,« K (Af re + 4P i>KK ) 



(19) 



(20) 



where Pi !fiv — (d^ipi\Vi\d v ipi), and the operator P s ; = I — 
\ipi){tpi\ projects into the space perpendicular to \tpi). 

Now consider the strategy in which measurements are made 
by projecting pg onto the eigenvectors of L K , thus saturating 
the quantum Cramer-Rao bound for estimation of 8 K . Gener- 
ally, however, the elements of the Fisher information matrix 
for estimation in this basis 7 M „ (6* K ) deviate from the quantum 
bound I^v, 

Ipv (Q K ) = ^Pi (df* logpA logp, + ^'^^'^ j ) 

(21) 

where the first term represents 'classical' information, and the 
second term makes the 'truly quantum' contribution. For com- 
parison the corresponding quantum Fisher information matrix 
is given by 



Pi (<9 M logp^ log pi + 43£P iiM1/ ) 



(22) 



Derivation of Fisher information matrices 

The state is given by Eqns.(17) and (18). The Fisher infor- 
mation matrices for states of this form are given by (21) and 
(22). In what follows, 



£,r,s = ^2 X t blt r , ^ r = ^2 Cr,s = ^ X t t r , 



(23) 



are the moments of the coefficients of Xt = |at| 2 , with the 
constraint ^™_ x t — 1, and 6* is the binomial factor defined 
in the last section. We begin by reexpressing the probabilities 

as 



Pi = Pi 



(in\i>i) = E Kl 2 b i ^J2 x " b i = (24) 



k=l k=l 

Evaluating the relevant inner products gives 

\ 2 hk 



; / |q M ■ k \ a k\ bf 



k=l 



Pi 



(25a) 



(d^tpild^ipi) 



"k 2 \a k \ 2 b k _ 6,i 



k=l 



PlO 



(25b) 



(Md v ^i) =^^e" ifc M -LB, a k e ik *J -L 

k=i V Pi y V P { 



(25c) 



This reflects the geometry of the estimation problem, wherein 
the i] derivative is orthogonal to the initial state, 



{dr,ipi\dr,ipi) =^2d rl 



f -ik<t> J b k\ / akG ik<t> h k 



k=l 



k=l 




bf \ &,t£a,i - £i , 



(25d) 



(d^ild^t) =-iY / ka* k e- ik *J b Ld v l ak e ik *J$- 



(25e) 



Writing these in the shorthand P k ,u,v = (d^ipk\d u ipk) ~ 
(d^ipkltpk) (ipk\dvipk) introduced in the last section, we find 



Pl.rhrh = 



&,i€o,i „ &,i£o,i - %i I 



£ 2 
so, 2 

&,i£,o,i - l 



j Pi, mi 



(26a) 



(26b) 



Also required is the Fisher information of the distribution 
{pi}, given by the terms (<9 M log pid u log pi). In order to eval- 
uate these one needs to use some properties of £ n j, including 



d v bf = bf 



I 



V »?(!-»?) 



1=0 

;2 ( 



1=0 



V 



12 c 

(1 — IT) 2 
1=0 K " 



— 



1—1] 



(27a) 
(27b) 
3 f+ i (27c) 

(27d) 



These relations can be used to evaluate 



n Tl 1 

(dr, logpi d v logpi) = V — {d n pif = V — (d v £ ,i) 
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The other products are much simpler and follow triv- 
ially from d^pi = 0. That is, (d<j, logpi \ogpi) = 
0, {dtj, log d n \ogpi) = 0. Although derived for the loss and 
the phase parameters here, this general feature for any pair of 
parameters corresponding to hamiltonian and dissipative dy- 
namics. 

Combining these results, the quantum Fisher information 
matrix elements from Eq. (22) are thus given by 

n 

^p i ((^logp i ) 2 +45RP / , 00 ) 

1=0 

n / n f 2 \ 

n 

= E Pi {{dnlogpi) 2 + m.P Um ) 



Asymptotic profile 

This section provides a step by step derivation of the 
asymptotic results in the high loss, high flux regime. We be- 
gin with the expressions for the quantum Fisher information 
for phase and loss estimation given by Eqns. (29a) and (29b). 
We first separate into its constituent terms after having 
shifted the mean by a constant 

n N / t \ 2 



Irhrh — 4 



L 1=0 k=l 



n 1 

£- 



I 



E & ? 

k=l 



1-T) 



\a k \ 



(34) 



1=0 



— if* 

r)(l-r)) t] 2 \ 2 fW ,z 



1 



E 

1=0 



By considering only the first term, and noting that the second 
term is a positive quantity, one obtains an upper bound for 



V 2 £o,i 



I 



1-77 



1—7? f— ' 



k=0 



77(1 - 77) 

Xfj(f) — X^yj 

When evaluating the Fisher information matrix elements 
given by Eq. (21) for the estimation using the optimal mea- 
surements for <fi, we obtain the same forms for 1 except I m 
where the useful cancelation which occurs in Eq. (29a) no 
longer occurs. 

*-<££<^-<w-«(*-£fO 

*(«^--i( S2 -gg) ,30.) 

(;w =/ >)-° < 30b ) 

Also of interest is estimation using the optimal measurements 
for 77. One finds that since Pi :V 4> purely imaginary, and pi 
carries no <$> dependence, so 1^1 = ifjP = 1^2 = 0, and 



(29b) 
(29c) 



EE*? 

7=0 k=l 

In the limit of large n, it is reasonable to approximate the 
sums by integrals with k = nz — n(l — y) and |afc| 2 = 
a 2 (z)dz = — f 2 (l — z)dz = f 2 (y)dy. This implies x k = 
\a k \ 2 ^^a 2 (k/n)^^f 2 (l-k/n) 

Efc|«fc| 2 « / k\a k \ 2 dk = n / (l-y)f 2 (y)dy. (36) 
fc=o J ° J ° 

Under this transformation the normalisation condition trans- 
forms as follows 

i = ^iv>) = £ki 2 « / \a k \ 2 dk= / 

o Jo Jo 



k=0 



\a k \*6k= / f{y)dy. 
10 Jo 

(37) 

Hence by Eqns. (35) and (36), the first term of Eq. (34) is 

n „i 



1-7? 



£&ki 



fc=0 



77 7/ 

1-vJo 
nq 



(l-y)f 2 (y)dy 



1 



/£, = 

w t?(1-»?)" 



(31) 



A final result that is of some significance is the condition for 
commutativity of SLDs [31]. In terms of the shorthand intro- 
duced in the last section, it is 

n 

Tt[p(L^- L U L^)} = AiJ^PftPl,^, (32) 

1=0 

which for the case of simultaneous estimation of 7/ and 0, us- 
ing the results of (24), (25) and (26), is 

n 

Tr [p (L v L,p - L^L n )] = 4 E (fc,l£o,l - 



I/^(v)dtf . (38) 
1 - ? / V Jo J 

Now consider the normalisation constants pi = 

Xm=Z l a fe| 2 - The continuous limit is taken using the 

de Moivre-Laplace theorem. Under the added condition that 

I R3 fc(l — 77) as fc — » 00 



6f 



1 



V2tt/c77(1 - 77) 
This leads to 

1 

Pi 



exp 



: exp 



;=o 
2 

-i — 
»7 



n el * 



1 ^2i:krj(l - 77) 

1 /' 



x n(l - r))V2n(T 2 V 2 



exp 



2krj{l — 77) 

{l-k{l- V )) 2 
2fcr?(l - 77) 

M - M) 2 



z 2ct 2 



(39) 



\a k \ dk 



a 2 (z) dz 



(33) 



1 



X n(l — rj)^/2na'- 



: exp 



(g ~ \tt 
2a 2 



a 2 (z) dz (40) 
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where /i = n ^_^ ■ & 2 — ^n^p an( ^ ^ ^ ast approxima- Using the approximation given in Eq. (39) to obtain a contin- 

tion is near equality for \i > a. A Taylor expansion to the uous limit form for the second term in Eq. (34) leads to 

first order a 2 (z) — (a(fi) + a'((J,)(z — fi) + 0(z - /i) 2 ) 2 = 

a 2 (/i) + 2a(fj,)a'((i)(z — fi) + 0(z — /i) 2 substitution into 

the integral, and taking the limit to infinity leads to the first 

order term being zero by parity. Due to the sharpness of the 

distribution, the limits can be taken to infinity. 

^ a 2 (aQ f°" 1 / (z-,i) 2 \ , a 2 (aQ 



k=l 



i V y/2irkr](l - rj) 



exp 



(I-fc(l-»/)) s 



2fc?7(l - 77) 



lajfel dfc 



i 1(1- ?7)v / 27R7 f V z 



exp 



M ~m) 2 
'z 2<r 2 



(z — /i)a 2 (z) dz 



i (1 - r))\/2ira 2 



exp 



(g - M) 2 
2a 2 



(z — /i)a 2 (z) dz, 



(41) 



where the same values of \i — an d cr 2 = ^j—p are used. When the limit is taken to infinity in this integral one loses 

the zeroth order term of the Taylor expansion by parity and retaining the first order term 



k=l V 

1 f 00 1 



l«kl 

exp 
1 



(z-m) 5 



1 - V J-oo V2TU7 2 

2a([x)a' (n) 

1 - V J-oo V27TC7 2 

2a(/i)a'(/i)er 2 



2a 2 



(z — /j)a 2 (z) dz 



: exp 



2a 2 



(z — [i) 2 dz 



1 — T) 



(42) 



r 



The second term in Eq. (34) is then 



1 

Pi 



Lfc=2 v 
n(l — 77) 



V 



\Otk\ 



a 2 (n) 

12 



2a(/j)a'(/i) cr 



I 



1—7] 

An 



n(l — r/) J 1 — X] \n 2 (l — rj) 



It, 



Taking the sum over I to the continuous limit, one obtains 

1 2 

2 



(=0 * Lfc=Z 
cN 



I 



4n 



n(l — 77) / 1 — ry \n 2 (l — 77) 



d/ 



4ti 2 



= 4 



2 „i 



2 , 1 



z 2 a /2 (z)dz 



?7 

1 — 77 



1-1 



(1 - vrrWv, 



(44) 



Putting together Eqns. (38) and (44), one obtains an asymp- 
totic form for Eq. (34) 

•I 



(43) 



- 4 



nr, 



I-77 



1 



2 fi 



yj (y)dy 



(l-y) 2 f 2 (y)dy 



(45) 
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Taking the upper bound of both integrals to infinity, the lower 
bound of the second integral go to zero, and using that the 
distribution is sharply peaked around y ps in the second 
integrand one obtains the final form 



with 7 3 = (2 — i) Tt Q 4 ? ?7 \2 ; which is equivalent to minimiz- 
ing the integrand when provided w > 1/2. In that case, the 
Euler-Lagrange method leads to the differential equation 



1 — 77 

Using similar techniques 



4 " f 2 (y) ] d» 



n 



Therefore, 



7j(l - 77) 



1- / yr(i/)d» 




(46) 



(47) 



yf(y) = 1 



A 2 f{y) 

dy 2 



(50) 



subject to the normalisation condition and the initial condition 
/'(0) = 0. The solution is an Airy function, in particular 



V 7 |Co|Ai 2 (Co) V7 



(51) 



7 '''' Arf 1 ^ ~ {I -if) 2 



f' 2 (y)dy 



(48) 



The final task is to maximize a convex combinations of the 
Fisher informations for phase and in Eqns. (46) and (48) with 
a weight < w < 1 assigned toZ^. The two quantities how- 
ever scale differently with n and 77. Consequently, a maximi- 
sation of their combination is unduly weighed by one of the 
quantities. To circumvent this, the Fisher information for the 
two parameters should be rescaled by their respective max- 
ima, given by = 477,77/ (1 — 77), and X,™ ax = 71/77(1 — 77). 
The final quantity to be maximized is given by 



w 



xmax \ I rmax 

<S><S> vv 

/>oo 

1-/ (yf 2 (y) + 7 3 f 2 (y))dy 

Jo 



(49) 



where £q = —1.01879 is the first zero of the first derivative if 
the Airy function, i.e., Ai'(£ ) = 0. Reverting to the quantities 
of original interest, 



xu ps ^ Ai 2 

n 7 |Co|Ai 2 (Co) 



7 777 



and 



Co 



(52) 



- (1 - 7lCo|)2^T x , Irm - ^w aX - 



(53) 



